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We review the use of Wilsonian renormalization group methods for quantum field theories at finite 
temperature. The implementations within both real and imaginary time formalism is carefully 
discussed. In particular, the question of gauge invariance is addressed in detail. A recently proposed 
real time thermal renormalization group is then used to derive the RG flows for Abelian Higgs models. 

Further applications are outlined. 


I. INTRODUCTION 

Most of the physically interesting questions in thermal 
field theory are outside the domain of validity of pertur¬ 
bation theory. This is true not only for static quantities 
like the magnetic mass or the equation of state for a 
quark-gluon plasma, but as well as for dynamical ones, 
like the plasmon damping rate close to the critical tem¬ 
perature. The reason for this break-down are IR diver¬ 
gences, which are difficult to control perturbatively. 

The use of reliable resummation procedures seems 
therefore mandatory. The Wilsonian or Exact Renor¬ 
malization Group is precisely a tool that allows for a 
systematic resummation beyond perturbation theory. 

All current implementations have in common, that 
they use the “known” physics in the UV as the starting 
point (see fig. 1). This region -the upper right corner- 
corresponds to the bare classical action of the T = 0 
theory. The goal is to find the corresponding effective 
action of the soft modes at non-vanishing temperature. 
This region -the lower left corner- describes the IR limit. 
The resummation problem is thus the question about how 
these two regions are related. 

For vanishing temperature, the flows towards the IR 
only integrate-out quantum fluctuations, and are de¬ 
picted by the flows along the vertical boundaries (for the 
4(i and 3d limits, respectively). For T 0, one might 
distinguish essentially three scenarios. 

1. The dimensional reduction approach aims at relat¬ 
ing the Ad, T = 0 parameters to those of an effective 3d 
theory at T = 0. This reduces the problem to a purely 
3d one, that is, to the problem of integrating out only 
quantum fluctuations in 3d. The temperature enters via 
the initial parameters of the effective 3d theory. 

2. Flow equations in the imaginary time formalism 
are used to directly relate the Ad couplings in the UV 
region at T = 0 with the renormalized ones at T yf 0. 
Both quantum and thermal fluctuations are integrated 
out, as the imaginary time formalism does not distinguish 
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between them. The Euclidean flow equation in Ad can 
directly be used, with the standard prescriptions of the 
Matsubara formalism. This scenario corresponds to the 
flow along the diagonal in fig. 1. 


3d-UV ’dimensional reduction’ 4d-UV 



FIG. 1. The qualitative difference of renormalization group 
flows for theories at Hnite temperature. The arrows indicate 
the flow towards the infrared. 

3. Flow equations in the real time formalism are used 
to relate the Ad renormalized couplings at T = 0 with 
the renormalized ones at T yf 0. A prerequisite of this 
approach is of course the knowledge of the renormalized 
Ad couplings in the first place. Contrary to the imagi¬ 
nary time approach, this one allows the investigation of 
non-static properties. The flow equation only integrates- 
out thermal fluctuations, that is, it describes how modes 
with momenta around k come into thermal equilibrium 
at temperature T. This RG flow corresponds to the flow 
along the base line in fig. 1. 

The paper is organized as follows: The sections II 
and III aim at giving an introduction to the key aspects 
of Wilsonian RGs. Sect. II reviews the Euclidean for¬ 
malism, and in particular the implementation at non¬ 
vanishing temperature and the inclusion of gauge fields. 
Sect. Ill is reserved for a recently proposed real time im¬ 
plementation. Sect. IV presents an application of the 
latter to the U(l)-Higgs model, while sect. V contains 
the discussion and an outlook. 
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II. WILSONIAN FLOW IN EUCLIDEAN TIME 

In this section we will outline the Exact Renormaliza¬ 
tion Group approach to quantum field theories in Eu¬ 
clidean space time In particular, we discuss the 

implementation for gauge theories ja^, and the appli¬ 
cation to thermal field theories ffCHiy in the imaginary 
time formalism. 


A. Coarse graining 

The main problem of perturbative methods for field 
theories at vanishing or finite temperature can be linked 
to the problematic IR behaviour of massless modes in less 
than four dimensions. It is therefore mandatory to find 
a regularization for them. Let us consider the case of a 
bosonic field </>. A particularly simple way of curing the 
possible IR singular behaviour of its perturbative propa¬ 
gator Pcj, consists in replacing it by a cut-off propagator 

P4,ek{pVk^) ■ ( 1 ) 

Here, we introduced a function 0^ which depends on 
a yet unspecified additional momentum scale k. 0^ is 
meant to be a (smeared-out) Heavyside step function: for 
large momenta p ^ k it goes to one (no regularization is 
needed), while it vanishes (at least with for p k. 
Eor any fc > 0, the above propagator remains IR finite 
and could safely be used within loop integrals. Finally, 
however, we are interested in the A: ^ 0 limit in which 
the regulator is removed. Thus we need to describe the 
fc-dependence of the theory, which brings into life the 
Exact Renormalization Group. Originally, it has been 
interpreted as a coarse-graining procedure for the fields, 
averaging them over some volume k~‘^. In this light, the 
IR limit corresponds to averaging fields over bigger and 
bigger volumes, i.e. to the limit k 0. 


B. The exact renormalization group 

A path integral implementation of these ideas goes 
back to 0,^ , where a regulator is used in order to distin¬ 
guish between hard (p^ > fc^) and soft modes (p^ < k'^V 
A slightly different point of view has been taken in g] 
(see also [Q), where a smooth cut-off has been employed. 
Following these lines, one obtains an effective theory for 
the soft modes only. The starting point is the functional 

expWk[J]= J Vfjjexp (^-Sk[4>\+ j 

( 2 ) 

Here, (p stands for all possible fields in the theory which 
we shall restrict to be bosons, for simplicity. (The ex¬ 
tension to fermions is straightforward.) J are the corre¬ 


sponding sources, and St = S+^kS contains the (gauge- 
fixed) classical action and a quadratic term AfcS'[^], 
given by 

^ y <I>*{-P) Rk{p) (p{p)- (3) 

It introduces a coarse-graining via the operator Rk(p)- 
Performing a Legendre transformation leads to the 
coarse-grained effective action T k [</>], 

TklcP] = -Wk[J] - AkSiP] + Tr Jcj>, p = (4) 

oJ 

where the trace Tr sums over momenta and indices. It 
is straightforward to obtain the flow equation for T^, 
w.r.t. t = Infc/A (with A some UV scale). The only ex¬ 
plicit fc-dependence in (||) stems from the regulator Rk, 
thus 


ftrfc = iTr|GfcM (5) 

with 

denoting the full (field-dependent, regularized) propaga¬ 
tor. For the time being, the regulator function is kept 
arbitrary. It can be chosen in such a way (see below), 
that 


lim Tk = S (7) 

k—t'OO 

lim Tk=T . (8) 

fc^O 

Therefore the flow equation connects the (gauge-fixed) 
classical action S with the full quantum effective action 
r. Solving the path integral (||) (for A: = 0) is therefore 
equivalent to solving (|) with the initial condition (g) 
given at some UV scale. One might read this approach 
as a path integral independent definition of a quantum 
field theory. 

Note that the flow equation (j^) is exact: no approxima¬ 
tions have been employed for its derivation. This means 
in particular that (|) is non-perturbative - it describes 
unambiguously the resummation of all (quantum and/or 
thermal) fluctuations. 

Solving the flow equation necessitates, however, trun¬ 
cations and approximations. One can easily recover the 
perturbative loop-expansion, or resummations thereof. 
However, (||) allows for more elaborate expansion 
schemes which are not confined to regions of small cou¬ 
pling constants. Commonly used is the derivative expan¬ 
sion, or an expansion in powers of the fields. Deriving (j^) 
w.r.t. the Helds gives then flow equations for the higher 
order vertices, which parametrize the effective action. 
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C. The regulator function 


Let us be slightly more explicit about the regulator 
function R^. We will impose the following constraints on 
the regnlator Rk such that IR Hniteness, (@) and (|^ are 
ensured: 


(i) Rk has a non-vanishing limit for —> 0, typically 

like a mass term Rk ^ k'^. 

(ii) Rk vanishes in the limit fe ^ 0, and for ^ k^. 

(hi) For k ^ oo (or k ^ A with A being some UV scale 
much larger than the relevant physical scales), Rk 
diverges like A^. 


Condition (i) reflects the IR hniteness of the propaga¬ 
tor at non-vanishing k even for vanishing momentum p 
(which is why the regulator has been introdnced in the 
hrst place). Condition (ii) ensures that any dependence 
on Rk drops out for fc —s- 0 (that is to say that F^^o 
reduces to the full quantum effective action F), and that 
large momenta modes are suppressed, (ie. integrated- 
out) . From condition (iiL we conclude that the saddle 
point approximation to (g) becomes exact for fc —> A, 
and Ffe^A reduces to the (gauge-hxed) classical action 
S. The regulator function is related to &k in (|^ as 


0fc 



Rk{p^) 

p2 + i?fc(p2)- 


(9) 


One easily verihes that any Rk with the properties (i)-(iii) 
yields a (smeared-out) Heavyside step function, when in¬ 
serted in (i- We also conclude from (||) that the oper¬ 
ator dtRk is a (coarse-grained) 5-function. This is con¬ 
sistent with the flow equation (™. At any fixed scale k, 
only loop-momenta p"^ around ^ can contribute to the 
change of Ffc. All other momenta are suppressed because 
of dtRk having support only in the vicinity of k^. This 
is the essence of a Wilsonian philosophy based on the 
integration over inhnitesimal momentum shells. 

Typical classes of smooth regulators used in the liter¬ 
ature are exponential ones with 


Rk{p'^) = 


or algebraic ones with 


exp(p2/fc2)’^ — 1 


Rk{p'^) =P^ 


2 \ ^ 


( 10 ) 


( 11 ) 


For n = 1, both classes have a mass-like limit for small 
momenta. The limit n —> oo corresponds to the sharp 
cut-off limit l]^ . 

Two comments are in order. The first one concerns the 
case of a mass-like regulator Rk = k^. This regulator is 
somewhat special. First of all, it is independent of mo¬ 
menta. The operator dtRk in (|^) is -for this particular 
regulator- neither peaked nor sufficiently suppressed for 


large momenta. Thns, all momenta p do contribute to 
Ffe at any k {i.e. the second part of (ii) is not fulfilled). 
Typically it is observed that the convergence properties of 
approximate solutions to the flow equation are worse for 
a mass-like regulator then for exponential ones. Further¬ 
more, the unsuppressed large momenta modes introduce 
additional UV divergences. Although analytical compu¬ 
tations simplify tremendously with this regulator, it has 
to be taken with care. 

The second comment concerns the (in-)dependence of 
physical quantities on the shape of the regulator function. 
Obviously, physical observables should not depend on the 
particular regulator chosen ||^, and will not depend on 
Rk, if the flow equation is integrated down to fc = 0. 
This is an immediate consequence of (^. However, any 
computation has to resort to some approximation (typi¬ 
cally only a finite number of operators are considered for 
an expansion of Ffc). Thus, approximations can intro¬ 
duce a spurious scheme dependence. Any approximation 
scheme, that yields large quantitative or even qualitative 
corrections is not acceptable and has to be discarded. 
Therefore it is mandatory to compute the scheme depen¬ 
dence within a given approximation. At the same time, 
this is an efhcient way to check the viability of a given 
Ansatz fl^ . 


D. Gauge invariant Green functions 

What can be done for gauge theories? The obvious 
problem is that the regulator term -being quadratic in 
the helds- is not gauge invariant, which raised some crit¬ 
icism about the present approach. The qnestion arises as 
to which extent gange theories can be handled, although 
the regulator term (^ seems not to be compatible with 
gauge invariance. This problem has been considered us¬ 
ing the background field method 0,^ , modihed Slavnov- 
Taylor or Ward Identities or (perturbative) fine tun¬ 
ing conditions |^. 

We will follow 1^ to illustrate the problem, and to 
show that gauge invariance of physical Green functions 
can indeed be maintained. Let us consider the example 
of an SU(iV) gauge theory in an axial gauge, with (/) = 
and the gauge fixing 

= ^Tr ^ (12) 

Here, ^ is the gauge fixing parameter (chosen to be mo¬ 
mentum independent) and is a fixed Lorentz vector. 
Axial gauges have the nice property that the ghost sec¬ 
tor decouples. The problem of Gribov copies is therefore 
absent, and the number of Feynman diagrams is signif¬ 
icantly reduced. Furthermore, the spurious singularities 
as encountered within a perturbative approach have been 
shown to be absent . Finally, the axial vector which 
appears in the gauge fixing condition, has a natural ex¬ 
planation within thermal field theories as the thermal 
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bath singles-out a rest frame characterized by a Lorentz 
vector. 

Let us now perform an infinitesimal gauge transforma¬ 
tion. This leaves the measure in (j^) invariant and leads 
to a functional identity, the so-called modified Ward Iden¬ 
tity (mWI) = 0, with 

Wt\A\ = 

-a J (13) 

and -I- . The mWI contains all terms 

of the standard Ward Identity (WI), i.e. the first line of 
Additionally it contains also a term proportional 
to Rk, which is a remnant of the coarse-graining. We 
observe that this term vanishes for any regulator in the 
limits k oo and fc —> 0. In particular, it vanishes iden¬ 
tically for a mass-like regulator Rk = thus reducing 
the mWI to the usual WI in these cases. The flow of 
yyfc[A] is obtained from ®) and (p^: 

a.ws = -iiY (g.^g4 X MIJ . (14) 

The flow ( p^ has a fixed point Wfe^] = 0. Thus, if 
rfe[A] solves ( p^ ) at some scale fco and evolves subse¬ 
quently according to (||), then Tk[A],k < fco, fulfills 
as well >Vfc[A] = 0. In particular, Tfe^o^] will obey 
the usual WI, which establishes gauge invariance for the 
physical Green functions. 

This approach has been used to prove that the 1- 
loop /3-function of SU(A^) gauge theory coupled to Nf 
fermions is indeed universal, independent of the choice 
for Rk or ^ (in d = 4) ||J^ . 

It is worth stressing that the mWI plays a double role 
when it comes to approximate solutions, i.e. for a trunca¬ 
tion of rfc[A]. First of all, the mWI can be implemented 
even in these cases. Perturbatively, this is well known, 
and sometimes denoted as the Quantum Action Princi¬ 
ple. A general procedure that allows to respect the mWI 
for numerical implementations even beyond perturbation 
theory can also be given (for the details, see @). At the 
same time, the mWI allows the control of the domain of 
validity for a given truncation |l^ . This error control is 
quite welcome also on a computational level as it avoids 
to go to the next order in the chosen expansion. 


E. Imaginary time formalism 

It is straightforward to upgrade the above approach to 
the case of non-vanishing temperature within the imag¬ 
inary time formalism [|l0| . The flow equation contains 
a loop integral over some momentum dependent func¬ 
tional. Therefore, the only changes concern the Tr in 
the flow equation, which becomes a sum over Matsubara 


frequencies, and the 0-component of the loop momenta, 
which is discretized 



po 2'KnT . 


(15) 


Note, however, that the flow equation now connects the 
UV parameters of the dd theory at T = 0 with the IR 
ones at T ^ 0. Thus, both quantum and thermal fluctua¬ 
tions do contribute to the flow equation. This procedure 
corresponds to integrating along the diagonal as depicted 
in fig. 1. It has been applied to phase transitions in scalar 
EMI and gauge field theories @0. 

At this point it is interesting to note the similarity be¬ 
tween exact flow equations and the approach advocated 
in [T^ to compute the non-perturbative pressure. Indeed, 
the method of |^] can be seen as a flow equation with 
a mass-like regulator Rk = fc^. But as we commented 
earlier, a mass term regulates only marginally. In order 
to avoid the additional UV problem for large momenta 
fluctuations, one should consider instead differences, like 
P[T] — Up]. This makes these divergences to cancel out 
and yields a well-defined flow equation. Even more inter¬ 
esting is the extension to gauge theories which has not 
been studied yet. For a mass-like regulator the second 
line in ( 0 ) vanishes identically. This corresponds to the 
statement that gauge invariance can be maintained for 
any fc (in this particular case), and is a special feature of 
the axial gauge fixing used. The above leads to the con¬ 
clusion that the generalization of 0 to gauge theories is 
most conveniently done within axial gauges |^. A more 
detailed account is given in [0 . 


III. WILSONIAN FLOW IN REAL TIME 

The philosophy of the previous section is most appro¬ 
priate for static situations, that is to say for equilibrium 
physics, and can be used to compute physical quantities 
at first order phase transitions {i.e. free energies, surface 
tension, latent heat), or at second order ones {i.e. criti¬ 
cal exponents, equations of state, amplitude ratios). Yet 
a number of interesting physical problems are related to 
non-static properties of quantum field theories, and the 
question raises whether this approach can be extended 
to space-time with Minkowskian signature. 

Recently, a strategy has been proposed for integrating- 
out the temperature fluctuations within a real-time for¬ 
mulation of thermal field theory 0]. The key idea is 
to introduce a thermal cut-off for the on-shell degrees of 
freedom. This philosophy has several advantages. By 
construction, it allows to study precisely the effects of 
thermal fluctuations only, and it is not restricted to static 
quantities [0. As thermal fluctuations are on-shell, it 
is straightforward to guarantee a gauge invariant imple¬ 
mentation of this coarse-graining even for intermediate 
coarse-graining scale fc [Q. However, no statements can 
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be made regarding quantum fluctuations. They have to 
be included from the onset in the initial condition. 


A. Real time formalism 

At finite temperature the fields (j) s-re defined on a 
contour C in the complex time plane [^3|-p5|. In the 
real time formulation this contour consists of a forward 
branch (the real axis in the complex time plane) and 
a backward branch (parallel to the forward branch, at 
arbitrary distance a < 1/T). It is convenient to intro¬ 
duce two separate set of fields (f>i (the original fields) 
and 4>2 (the so-called thermal ghosts), living on the two 
branches. Thus there is a doubling of the degrees of free¬ 
dom and the propagators become 2x2 matrices, deter¬ 
mined by the boundary conditions. Let us consider the 
case of a scalar field of mass m. Its free propagator 
at r = 0 is 


— m? -|- ie 

The corresponding tree level real time propagator in mo¬ 
mentum space is 

= Q[D^{p)] + A^p) Nt{\po\) B. (17) 

We introduced A'^(p) = D'^{p) — (^D'^’{p)y and the ma¬ 
trices Bij = 1, {i,j = 1, 2) and 

/ A^9{-po)\ 

Q[Dy = (18) 

VA^0(po) / 

The function Nt denotes some thermal distribution func¬ 
tion, which, at thermal equilibrium, would be the Bose- 
Einstein distribution. It contains basically all the ther¬ 
mal information. In the e: ^ 0 limit, we observe 

A'^(p) —> —2in6{p'^ — n?) (19) 

and conclude that the thermal part of the tree level prop¬ 
agator involves on-shell degrees of freedom only. 



B. Thermal coarse graining 

The question is how the thermal propagator could be 
modified in order to implement a thermal coarse grain¬ 
ing. The proposal of consists in a scale dependent 
modification of the thermal distribution function, 

AT.fe = ATefc(IPlA) (20) 

The corresponding tree level cut off thermal propagator 
obtains from (0) through the replacement (p0|), i.e. 

D^k = QlD^^ip)] + A^{p) iVT.fcdpol, IpI) B. (21) 


Eq. (|2^) may be interpreted as the thermal analogue of 
(|^. The thermal distribution Nt is switched on mode 
by mode through the ©^-function, whereas in (H), the 
propagation of longer wave length modes is switched on. 
The modes with |p| 3> fc will be in thermal equilibrium 
at temperature T, while those with k 3> |p| remain in 
equilibrium at the temperature T = 0. We shall use in 
the sequel a sharp cut-off 

0fe(|p|)=0(|p|-fc). (22) 

Alternatively, an exponentially smooth cut-off is given by 

glPol/T 

~ ^ 1 _l_ g|p|/fe(gbol/'r _ 1) ■ 

Note that (^ ) yields a modified Bose-Einstein distribu¬ 
tion very similar to the one proposed by Nair -although 
within a different prospective- in . 


C. The thermal renormalization group 

A path integral formulation of these ideas has been 
given, following in full analogy to the Euclidean case 
pOf . Let us first introduce sources Ji {i = 1,2) for the 
fields (pi in order to obtain the path integral representa¬ 
tion of the generating functional of real-time cutoff Green 
functions as 

Zk[J] = y expi(iTr(/)i ..'/’I 

+5i„tM+TrJ,</),). (24) 

The trace corresponds to the sum over all fields and in¬ 
dices (only the thermal one have been given explicitly), 
and momentum integration. 5'int[(/)] is the bare interac¬ 
tion action 


iS'int)*/’] — ‘5'int['A2]. (25) 


The flow equation for Zk[J] can be derived easily from 
(p^) and reads, using t = ln(A:/A), 


dtZk[J] = --Tr 




(26) 


We define as usual the cutoff effective action Lfc as the 
Legendre transform of the generating functional of the 
connected Green functions, using IEfc[J] = ilnZfc[J], as 

Lfc [(P] = Tr J,cP, - i p, ^3 - Wk [ J] (27) 

with 


pi 


6Wk[J] 

SJ^ 


(28) 


where we have isolated the free part of the cutoff effective 
action and used for the classical fields the same notation 
as for the quantum fields. 
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The flow equation for the cutoff effective action 
follows as 


Tfc [(/> = const] 



(31) 


dtrkm = ^Tr 


dtD 


-1 

T,fc 



6(j)S(l) ) 


(29) 


(thermal indices have been suppressed now). This flow 
equation is the thermal analogue of i). Given that the 
initial condition for Ta [^] at A T is the full renormal¬ 
ized theory at zero temperature, the above flow equation 
describes the effect of the inclusion of thermal fluctua¬ 
tions at a momentum scale around jpj = k. At any fixed 
scale fc, Tfc describes a system in which only the high fre¬ 
quency modes jpj > /c are in thermal equilibrium, while 
the low frequency modes jpj < k do not feel the thermal 
bath and behave like zero temperature modes. 

The flow equations for all the possible vertices are ob¬ 
tained by expanding ( p^ ) in powers of the fields. For 
their derivation it is helpful to note that the cutoff effec¬ 
tive action has a discrete Z 2 symmetry 

Tk[cPuh] = -nm,ri], ( 30 ) 


We shall evaluate it for a field configuration (p = 4^ which 
gives a non-vanishing v.e.v. ^ only to the real part of one 
component of (j)i. All other fields are then set equal to 
zero. In |2^] it was shown, that the Euclidean potential 
(for one real scalar field) is given by 

dVki^) ^ dVkjcP] 
dcj) d(j>i 

This relation can be generalized to an arbitrary number 
of fields 1^]. Apart from an irrelevant constant, the ef¬ 
fective potential is 

Vk[$) = (</.), (33) 


(32) 


where 





(34) 


which relates different vertices to each other. 


IV. APPLICATION: U(l) HIGGS THEORY 


denotes the tadpole. Thus the flow equation for the Eu¬ 
clidean potential is deduced from the one for the tadpole, 
using (^. 


In this section we will apply the real-time thermal RG 
to an Abelian Higgs model. A derivation of the flow 
for the effective potential is given. A more elaborate 
presentation will be given elsewhere . 

The U(I) Higgs model is quite appealing as a testing 
ground for the feasibility of this approach. Eirst of all, it 
is an important model for cosmological phase transitions. 
Furthermore, the starting point for the use of (p^, that 
is the renormalized action at vanishing temperature, is 
well known and can be computed with high accuracy. 
This model has several different mass scales, which we 
expect to thermalize and decouple at different coarse- 
graining scales. Finally, in the limit T —> 00 we expect 
to hnd flow equations for the purely id Abelian Higgs 
model. This might shed new light on the superconducting 
phase transition in id. 

A. The Euclidean effective potential 

The field content of this model is given by N complex 
scalar fields, an Abelian gauge field, ghost fields, and 
their thermal partners. We will use the Landau gauge 
throughout. For the time being, we shall be interested 
in the phase transition at finite temperature, whose de¬ 
tails are encoded in the coarse grained effective potential. 
Therefore, we have to relate the Euclidean effective po¬ 
tential Vk{4>) to 14 [^]. Evaluating the cutoff effective 
action for constant field s results in 


B. Approximate flow equations 

In order to obtain a closed set of flow equations, we 
have to employ some approximations regarding higher 
order vertices. We shall employ the leading order ap¬ 
proximation in a derivative expansion. This implies that 
the wave function renormalization of the scalar field is 
neglected. Eurthermore, we shall assume that the action 
remains at most quadratic in the Abelian gauge field. 
The (possible) imaginary parts of the scalar and photon 
self energies are neglected. 

These approximations allow a closed set of flow equa¬ 
tions for the functions 14(p) and Uk{p), where p = (jxp*, 
and Uk{p) is the (field dependent) coefficient in front 
of the operator in the action. The flow for Uk is 
related to the longitudinal part of the photon self en¬ 
ergy at vanishing momenta. For the time being, we dis¬ 
card the distinction between the temporal and the spa¬ 
tial gauge field component. (The more general Ansatz 
C/fcA2 ^ C/i.fcAg + U2Mf is able to correctly describe 
the decoupling Debye mode.) It is useful to introduce 
the following shorthand notations 


Di(p) = + 2pV''(p) 

(35) 

n2(p) = ^k^ + V'(p) 

(36) 

^3(p) = v'fc2 + 17fe(p) 

(37) 


in terms of which the flow for 14 (p) reads 
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dtVkjp) 

Tk^ 


- ^ (-^s/r)] 6{nl) 

- Ain[l-exp(-03/r)] 0(0^). (38) 

An analogous flow equation is obtained for the function 
Uk{p). Now we will rescale all the dimensionful quanti¬ 
ties with appropriate powers of T and k in order to ob¬ 
tain dimensionless flow equations for the iV-component 
Abelian Higgs model. To that end, we shall introduce 
the following functions (using ' = dp): 


v(.p,t) = V{p,t)/Tk^ 

(39) 

u{p,t) = U{p,t)/k'^ 

(40) 

p = p/kT, 

(41) 

n,/k (^ = l,2,3) 

Wi {p,t) = \/l + v' + 2pv" 

(42) 

W 2 ip,t) = \/l + v' 

(43) 

W3{p,t) = Vl + u. 

(44) 


Within this notation the flow equation for v takes the 
form 


dtv = —3v + pv' 

In [1 - exp {-wik/T)] 

- 1) In [1 - exp {-W 2 k/T)] 

-3^^^1n[l - exp(-W3fc/r)] (45) 

Note that the flow still explicitly depends on the ratio 
k/T, which is just a consequence of the presence of two 
independent scales. 


C. Low and high temperature limits 


dtv = 






dTT^ 

(2iV - 1) ln[l-t u'] - 3 


2pv"] 


ln[l 


(46) 


Note that we have suppressed terms proportional to 
9{'uP‘)h\{k/T)^ as their contribution for u>^ > 0 is field- 
independent. 

It is interesting to compare (^^ with the flow equation 
directly derived in 3d via the effective average action ap¬ 
proach [|| . If we specify a sharp cut-off regulator and ne¬ 
glect the scalar anomalous dimension, the corresponding 
flow equation for the Euclidean dimensionless potential 
reads ll^,0 


dtv^ = -3v^ + “ 4 );^ 


■^Pv'^ 
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-^^Wl+4]-^ln[l + 2eV] 


(47) 


where jk denotes the dimensionless gauge cou¬ 

pling squared. 

Two comments are in order. First of all, the role of the 
(field-dependent) mass of the gauge field 
has been taken over by the function u, with = 
u{p)k‘^ = Uk{p)- With the computation of Uk{p) we 
would obtain therefore the full field dependence of the 
Abelian charge. In it was argued that the field de¬ 
pendence of might play an important role close to the 
critical points. Note also that the numerical coefficient 
in front of the last term in (|4^) differs from that in (H) 
(he. 3 instead of 2). As mentioned earlier, this is due to 
the Debye mode, which, in the present approximation, 
can not decouple properly, and therefore still contributes 
to the flow in the high temperature limit. 

The second point concerns the 0-functions, absent in 
(^). They ensure that the flow (46) will not run into 
a singularity nor develop an imaginary part. As soon as 
the arguments in the logarithms tend to negative values, 
the 0-function cuts their contribution off. In (^), this 
is not so obvious. While solving (^), however, one ob¬ 
serves that the flow does indeed avoid the singularities 
automatically [p9|-^. 


By construction, the flow equations (for v and u) only 
integrate-out the temperature fluctuations. This is why 
the initial conditions -given for the functions u and v at 
some large scales k- already have to be the full quantum 
effective parameters of the 4d theory at vanishing tem¬ 
perature. Therefore, the flow equations have to vanish in 
the low temperature limit since no further fluctuations 
need to be taken into account. That this is actually the 
case is easily seen in (^5|). For T ^ 0, the exponential 
factors exp(—wfc/r) k/T suppress any non-trivial flow. 

The high temperature limit is much more interesting. 
For T ^ oo we would expect to recover the purely three- 
dimensional running of the couplings. Expanding ln[l — 
exp(r(; k/T)] = In w -f ln(fc/r) -|- 0{w k/T) we obtain 


V. DISCUSSION AND OUTLOOK 

We reviewed main features of the Wilsonian RG and 
argued that they represent a systematic and efficient 
tool for applications to thermal field theories. In par¬ 
ticular, we emphasized that even gauge theories can be 
handled systematically in both the real and imaginary 
time formalism. The particular differences between these 
two approaches have been discussed, the main one be¬ 
ing that the imaginary time formalism is adequate for 
the computation of static quantities, while the real time 
approach allows the study of non-static quantities and 
non-equilibrium situations. 
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The application to an Abelian Higgs model extends 
previous applications to an interesting toy model for cos¬ 
mological phase transition. A detailed study of the hrst 
and second order phase transition is now feasible. This 
will also allow for an independent check of the pertur¬ 
bative dimensional reduction scenario, which is at the 
heart of recent Monte Carlo simulations. The extension 
to the electroweak phase transition seems to be straight¬ 
forward, although more elaborate approximations have 
to be employed in this case. It would be particularly in¬ 
teresting to study the critical points of both models. In 
the U(l) case, one expects two critical points (describing 
the second order phase transition, and the end point of 
the first order phase transition region). This approach 
might even open a door to a better understanding of the 
superconducting phase transition, which corresponds to 
the large T limit. In the SU(2) case one expects to find 
only one critical point, the end point of the line of first 
order phase transitions. This hxed point was recently 
discovered to belong to the Ising-type universality class 
A held theoretical determination of the end point 
and the corresponding critical exponents is still missing. 
Up to now only Monte Carlo simulations have been able 
to study this parameter range. 

With these tools at hand a number of other interesting 
problems can now be envisaged. An open question con¬ 
cerns for example the thermal /3-function of QCD, which 
has been computed by a number of groups, with remark¬ 
ably different results (see and references therein). 
The Wilsonian RG, and in particular the heat kernel 
methods used in 0- can be employed even within the 
imaginary time formalism and should be able to resolve 
this point. It seems also be possible to construct a gauge 
invariant thermal renormalization group within real and 
imaginary time along the lines indicated earlier . This 
would be a very useful extension of |l^] to fermions and 
gauge helds. 
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